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A model for the formation of supermassive primordial black holes in galactic nuclei with the 
simultaneous suppression of the formation of intermediate-mass black holes is presented. 



^ ■ A bimodal mass function for black holes formed through phase transitions in a model 

with a "Mexican hat" potential has been found. The classical motion of the phase of a 

complex scalar field during inflation has been taken into account. Possible observational 

O ■ manifestations of primordial black holes in galaxies and constraints on their number are 

. discussed. 
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£ ■ INTRODUCTION 



There is almost no doubt that the nuclei of most structured galaxies host supermas- 
sive black holes (SMBHs) with masses ~ 10 6 - 1O 1O M (Antonucci 1993; Falcke and Hehl 
2003; Ferrarese and Ford 2004). There is also observational evidence for the presence of 

CO 

in 



intermediate-mass (~ 10 3 — 1O 5 M ) black holes (IMBHs) in galactic disks and globular 
star clusters (Gebhardt et al. 2002; Lanzoni et al. 2007; Noyola et al. 2008; Strohmayer 
and Mushotzky 2009). In contrast to stellar-mass black holes produced by the explosions 
of massive supernovae, the origin of SMBHs and IMBHs remains enigmatic. Several sce- 
narios for the formation and subsequent growth of SMBHs have been proposed (for reviews 
and references, see Rees 1984, Begelman et al. 1984, Dokuchaev 1991, and Dokuchaev et 
al. 2007), which can arbitrarily divided into two groups: astrophysical and cosmological. 
In astrophysical scenarios, IMBHs and SMBHs result from the dynamical evolution and 
gravitational collapses of gas clouds or of central star clusters in galactic nuclei or from 
multiple mergers of stellar black holes. Cosmological scenarios suggest a very early for- 
mation of primordial black holes (PBHs) at the radiation-dominated stage from adiabatic 
density perturbations or through phase transitions in the early Universe. A cosmological 
scenario for the formation of massive PBHs through the collapse of isothermal baryonic 
charge fluctuations that arise during phase transitions at the inflationary stage has also 
been proposed (Dolgov and Silk 1993; Dolgov et al. 2009). 

The currently most popular models are those of hierarchical merging where SMBHs 
are formed by multiple mergers of less massive IMBHs (Miller and Colbert 2004; Kawakatu 
et al. 2005). The latter could be formed, for example, during the collapses of gas clouds in 
small dark matter halos at pregalactic epochs, during hierarchical merging of stellar-mass 
black holes that are the remnants of the first stars (Islam et al. 2003), or during accretional 
swallowing of the first massive stars by PBHs in their interiors (Bambi et al. 2008). 
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One of the scenarios for the formation of PBHs through phase transitions was de- 
veloped by Rubin et al. (2000, 2001a, 2001b), Khlopov and Rubin (2004), Khlopov et 
al. (2005), and Dokuchaev et al. (2005, 2008). Apart from the presence of SMBHs in 
galactic nuclei, this scenario also predicts the existence of IMBHs in galactic halos at a 
considerable distance from their centers. The distinctive features of the scenario are the 
formation of black holes already at the radiation-dominated stage and a characteristic spa- 
tial distribution of black holes in clusters. Whereas "astrophysical" black holes are formed 
in potential wells at the centers of dark matter halos (Dokuchaev et al. 2007), PBHs 
are capable of producing an induced dark matter halo around themselves that appears as 
a dwarf galaxy (Dokuchaev and Eroshenko 2001, 2003; Dokuchaev et al. 2008) with a 
sharp increase in density toward the center. The mergers of IMBHs formed through the 
mechanism proposed by Rubin et al. (2000, 2001a, 2001b), Khlopov and Rubin (2004), 
and Khlopov (2005) and the formation of a peculiar class of galaxies and quasars around 
primordial SMBHs were considered by Dokuchaev et al. (2005, 2008). 

It is hard to imagine that absolutely all IMBHs will merge and produce central 
SMBHs in presentday galaxies. It is more likely that only some fraction of them, ~ 1, will 
merge, while the remaining IMBHs form a population of IMBHs in galactic halos. A similar 
situation arises in many models with primordial SMBHs. Indeed, if some mechanism 
predicts the formation of one primordial SMBH in the bulk of one galaxy, then it is hard 
to avoid the presence of a large number of less massive IMBHs in the same galaxy. In 
addition, even if SMBHs in galactic nuclei are not related by their origin to IMBHs, the 
IMBH mass functions tend to increase toward low masses. 

Thus, a common feature of the above SMBH formation mechanisms is the predicted 
abundance of IMBHs. Nevertheless, observational data constrain significantly the possible 
number of IMBHs in galaxies. Only a small number of objects can be attributed to 
IMBHswith some probability, for example, IMBHs in some globular star clusters. However, 
as yet there is no unequivocal proof of the presence of IMBHs at their centers, because the 
available observational evidence for their existence is questioned and comes under justified 
criticism (Illingworth and King 1977; McNamara et al. 2003; Baumgardt et al. 2003; 
van der Marel and Anderson 2009). Besides, the formation of IMBHs in globular clusters 
is severely hampered by the heating of their central parts by binary stars (hard pairs) 
formed during dissipative two-body encounters of single stars (Ozernoy and Dokuchaev 
1982; Dokuchaev and Ozernoy 1982). 

IMBHs could accrete gas and manifest themselves as ultraluminous X-ray sources 
(ULXs), those with a luminosity Lx > 10 39 erg s -1 . However, various models have been 
proposed to explain ULXs observed in the disks of some galaxies even without invoking 
IMBHs (Roberts 2007; King 2008). We have to say that the currently available observa- 
tional data give no unambiguous answer to the question about the existence and abundance 
of IMBHs in the Universe. There are only specific examples and some evidence and the- 
oretical arguments for their presence. In this paper, we propose a model in which the 
IMBH formation is suppressed under effective SMBH formation in galactic nuclei. This 
scenario can become topical if future observations will constrain significantly the number 
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of IMBHs. We propose a modification of the mechanism initially suggested by Rubin et 
al. (200 fb) that allows not only the presence of SMBHs at galactic centers but also the 
suppression of IMBH formation to be explained. 

THE FORMATION OF PRIMORDIAL BLACK HOLES WITH 
ALLOWANCE MADE FOR THE CLASSICAL MOTION OF THE PHASE 

OF A COMPLEX SCALAR FIELD 

The subsequent discussion is based on the primordial SMBH formation model de- 
veloped by Rubin et al. (2000, 2001a, 2001b), Khlopov and Rubin (2004), Khlopov et al. 
(2005), and Dokuchaev et al. (2005, 2008), who showed that if the inflaton potential has 
a local maximum, then the SMBH formation probability at the post-inflationary epoch is 
high. These authors described in detail the formation of a black hole through quantum 
field fluctuations near this maximum. The consideration is based on the shape of a po- 
tential known as a "tilted Mexican hat" , where the local maximum or, more precisely, the 
saddle point is located at points 9 = n, 2tt, ... (9 is the phase of a complex scalar field). This 
potential, which is used, for example, to describe inflation on a pseudo- Nambu-Goldstone 
field (Dolgov and Freese 1995), is 

V(|$|) = A[<r$-/ 2 /2] 2 , (1) 

with the minimum at ($) = fe^^/\/2. However, we will emphasize that the field $ under 
consideration is additional to the inflaton one. 

The problem is slightly simplified on energy scales smaller than /, because the mas- 
sive radial mode of $ (m ra d = A 1 / 2 /) can be discarded. The remaining light degree of free- 
dom is defined by the angular variable 9 = 0/ / that acts as a pseudo-Nambu-Goldstone 
boson (PNGB). Through spontaneous U(l) symmetry breaking on energy scales ~ A <C /, 
the PNGB potential acquires an additional term, 

V{9) = A 4 [l -cos0], (2) 

In the equation of motion for the field 9 with potential (1), (2), 

• A 4 

9 + 3H9 + — sin# = (3) 
/ 

where the Hubble parameter H remains constant during inflation, because the field 9 is 
not the inflaton one. During inflation, the quantity 9 changes slowly and only after the 
end of the inflationary stage does 9 execute rapid oscillations near the minima necessary 
for efficient production of matter particles and heating of the Universe. Depending on 
initial conditions, the field 9 can roll down to one of the minima, 9 min = or 9 min = 2n. 

During the inflationary stage, when the field 9 is in the valley of the "Mexican hat" 
|$| = f/y/2, the space is divided into many causally disconnected regions. The values 
of 9 in these regions slightly differ due to quantum fluctuations. In some of them, the 
field 9 may turn out to be on the other side of the maximum 9 max = ir and after the end 
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of inflation may roll down to a value of the minimum different from that to which most 
other regions tend. Thus, a Universe that consists of chaotically distributed domains with 
fields of or 2ix inside emerges after inflation. The neighboring domains are separated by 
the field walls whose subsequent evolution leads to the formation of PBHs. The number 
of formed PBHs and their mass depend strongly on the tilt of the potential A and the 
symmetry breaking scale / at the onset of inflation. Dokuchaev et al. (2005) chose model 
parameters at which clusters consisting of initially massive PBHs, the largest of which 
reached a mass ~ 4 • 1O 7 M , were formed. Such massive PBHs can subsequently serve as 
the nuclei of protogalaxies, increasing their mass to ~ 1O 9 M through accretion. However, 
in addition to massive PBHs, a large number of intermediate-mass (~ 10 4 M Q — 10 6 M Q ) 
PBHs are formed. As was discussed in the Introduction, their number may turn out to 
be too large, which will be in conflict with observational data. Themechanism for the 
suppression of IMBH formation considered here is based on the following almost obvious 
fact. The farther the initial phase at which the present-day Universe was formed from 
the maximum (in our case, 9 — ir), the smaller the volume of space will be filled with a 
phase > 7i. Indeed, after each e-fold, spatial regions with phases advanced toward the 
maximum of the potential by some small 89 appear due to field fluctuations. The time, 
the number of e folds needed for the maximum of the potential to be reached, increases 
sharply with increasing difference between the initial phase and 9 = ir. Therefore, the total 
number of closed walls and, hence, black holes formed from them also decreases sharply. 
In other words, the dispersion of the phase distribution increases with time and the farther 
the mean field from 9 = ir, the fewer the black holes will appear subsequently. If the mean 
phase recedes from 9 = tt with time, then the formation of IMBHs will be suppressed. 
This is possible if the classical motion of the phase is taken into account. 

So far the classical motion of the phase 9 of a complex scalar field at the stage of 
inflation has been neglected. This is true if the tilt of the potential V(9) defined by the 
parameter A is negligible. If, however, this assumption does not hold, then, as a result of 
the classical motion, the mean field recedes from the maximum and the IMBH formation 
probability turns out to be suppressed. Since themain idea of our paper is take this effect 
into account, let us consider it in more detail. We will measure the time in e-folds N 
before the end of inflation. Let there be some region of space appeared N r e-folds before 
the end of inflation and filled with a field with phase 9 r (N r ). 

We will be interested in the phase change with time through both quantum fluctu- 
ations and classical motion. An appropriate mathematical apparatus was developed by 
Rey (1987), who proposed representing the field as the sum of the classical term and 
fluctuations $ about it, 

9 = e(t)+&, tf<0 (4) 

The probability density P/(i?) for detecting the fluctuational part of the phase d satis- 
fies the Fokker-Planck equation, whose solution has the form of a Gaussian distribution 
(Starobinsky 1982). Obviously, the phase distribution of interest to us can be obtained by 
the simple substitution P{9) = P f (Q - 9). 

The tilt of potential (1), (2) is small at the chosen parameters A = 1.75H and / = 
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10H and we will take it into account only when calculating the main, classical contribution 
0(t) in (4). The small corrections due to the fluctuations d in (4) will be taken into 
account by neglecting the tilt of the potential lest the accuracy of the calculations be 
exceeded. In addition, analytical results are known for the spatial phase distribution (see, 
e.g., Starobinsky, 1982; Rey 1987). 

In terms of e-folds, the probability of finding the phase in the interval (9, 59) at a 
certain point of space is (Khlopov and Rubin 2004; Rubin et al. 2001b) 



SP(9) = 



1 



2% (N r - N) 



exp 



(9 r (N) - 9f 
259 2 (N r — N) 



59 = 



H 



(5) 



Here, as was said above, we neglected the tilt of the potential. The classical part of the 
phase 6(t) is represented as 9 r (N), where the number of e-folds is N = Ht. 



The condition for the appearance of one, and only one, SMBH at time N 

5P(n) 



e 3(N r -N ) 



(6) 



indicates that black holes are formed when the extremum of the potential 9 = tt is crossed. 
Now, from the equation 



3 3(Af r -JV ) 



= J2ir (N r - N ) exp 



(6r(N ) 



IT 



259 2 (N r - JV ) 



(7) 



we will find the phase 9 r (N ) from which the spatial region must be produced for the 
SMBH formation conditions to emerge N e-folds before the end of inflation, 



9 r (N ) = 7T - 59J(N r - N ) [6 (iV r - iVo) - ln27r (iV r - No)]. 



(8) 



Let us now determine the formation probability of lower-mass black holes produced 
at an e fold N, so that N < N < N r . Figure 1 show the dependence of the phase 9 on 
time expressed in e-folds obtained by numerically solving Eq. (3) and in the slow-rolling 
approximation. The excellent agreement indicates that the system is in the regime of slow 
rolling. In the subsequent calculation of the probability of quantum fluctuations in field 
9, it is appropriate to use not the numerical solution of Eq. (3) but its linear fit. Since 
intermediate-mass PBHs are formed from large domains emerging in a small number of 
efolds, the solution of (3) is fitted in the initial segment (N < 20). It is convenient to fit 
the dependence by a linear function, 



9 r (N) ~ 9 r (N ) - a ■ (N - N ) 



(9) 



The probability for the emergence of IMBH formation conditions during an e-fold number 
N is then 



1 



y/2n(N r -N) 



exp 



(e r (N ) - a .(N-N Q )-n 
259 2 (N r - N) 



21 



(10) 
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Given the number of causally disconnected regions exp[3(N — N r )] formed since the for- 
mation of the spatial region N r , we can find the total number of IMBHs in this region. A 
detailed description was given by Khlopov and Rubin (2004), who, in contrast to our case, 
disregarded the classical motion of the field, i.e., a — 0. 

Figure 2 presents the PBH mass distribution for A = 1.75H and / = 10H. As we see, 
the IMBH formation actually turns out to be suppressed. The mass function is bimodal 
in shape: the region of low-mass black holes is separated by a dip from SMBHs. 

ASTROPHYSICAL IMPLICATIONS 

In the preceding section, we showed that the PBH mass spectrum could have two iso- 
lated regions. The maximum at the highest masses gives SMBHs that can become SMBHs 
in galactic nuclei, possibly, by additionally increasing their mass through gas accretion. 
Less massive IMBHs form a population of black holes in galactic halos. Since their for- 
mation is suppressed, these IMBHs are very rare and they can be virtually unobservable. 
Let us discuss in more detail the paths of black-hole evolution from these two sets. 

Under the influence of dynamical friction, a black hole that was initially in the halo 
loses its angular momentum and approaches the galactic center. However, the dynamical 
friction mechanism is efficient only for fairly massive objects, more specifically, objects 
with a mass of more than < 10 7 M Q sink to the galactic center in the Hubble time. An 
important factor that contributes to the motion of black holes to the galactic center is 
the formation of an induced dark matter halo around the black hole (Dokuchaev and 
Eroshenko 2001a, 2003; Mack et al. 2007). The mass of the induced halo exceeds the 
black-hole mass approximately by two orders of magnitude. Thus, black holes with a 
mass > 1O 5 M surrounded by the induced halo accumulate and merge at the galactic 
center. If the internal structure of the induced halo depends only on the PBH mass, 
then the boundary of the induced halo is determined by the environment in which the 
halo is formed. The neighboring density fluctuations at a pregalactic stage will be an 
alternative center of attraction for dark matter starting from some distance from the black 
hole (Dokuchaev and Eroshenko 2001a, 2003). 

If a pair of IMBHs was formed at the galactic center and if a third IMBH approached 
the pair under the action of dynamical friction, then the slingshot effect is possible: the pair 
becomes even closer as energy is exchanged during the gravitational interaction, while the 
third IMBHacquires a high velocity and is ejected away from the galactic center (Haehnelt 
and Kauffmann 2002). The detection of such ejected IMBHs would serve as evidence for 
the formation of SMBHs from merging lower-mass IMBHs. Another method for testing 
various SMBH formation models is to search for gravitational wave bursts from black- 
hole mergers. For example, in the model from Dokuchaev et al. (2008), IMBHs in central 
clusters in galactic nuclei collide with one another with the generation of gravitational wave 
bursts (Dokuchaev et al. 2009), with the burst redshift distribution differing from that 
provided by "astrophysical" scenarios for the origin of SMBHs. Therefore, observations 
with the LISA telescope planned to be launched can provide a choice between various 
scenarios for the origin of SMBHs (Dokuchaev et al. 2009). 
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In an alternative scenario, for example, that considered here, the central black hole 
initially had a high mass. In the hierarchical picture of galaxy formation, a SMBH fell 
into a large galaxy together with one of the small protogalaxies and rapidly sank to its 
dynamical center through dynamical friction. The advantage of the model with massive 
primordial SMBHs is that it can successfully explain the presence of quasars at high 
redshifts (Dokuchaev et al. 2005). 

Lower-mass (< 1O 5 M ) black holes together with the surrounding induced dark 
matter halos have no time to sink to the center and remain in the galactic halo. They can 
be observable as ultraluminous Xray sources when passing through dense molecular clouds 
in the disk. Since the gas density in the galactic halo is low, accretion directly from the 
halo will not give a significant X-ray flux. However, the case where an IMBH moving in the 
halo carries a small induced dark matter halo and, most importantly, a cloud of baryons 
and an accretion disk is possible. In this case, the X-ray emission from the IMBH will 
be determined by the accretion of baryons from the induced halo. This scenario appears 
problematic, because no cooling of baryons, in particular, no star formation, takes place in 
low-mass halos. For this reason, there is no flow of gas onto the central IMBH. This also 
implies that the observed ultraluminous Xray sources, candidates for the role of IMBHs, 
are present only in galactic disks. 

Let us discuss the growth of IMBH masses through accretion in the time of their 
passages through the galactic disks. For Bondi-Hoyle accretion, a black hole has a lu- 
minosity L x = A-Krjc 2 G 2 M 2 pv- ? ' (Mii and Totani 2005; Mapelli et al. 2006), where v is 
the velocity with which the black hole moves through a gas with density p and rj is the 
mass-to-energy conversion efficiency This parameter is rj ~ 0.06 for disk accretion onto a 
Schwarzschild black hole and reaches 77 ~ 0.42 accretion onto an extreme Kerr black hole. 
Typically, rj ~ 0.1. In this case, the black hole mass increases with the characteristic time 
^bh = (4:irG 2 Mpv ~ 3 ) -1 . If the luminosity is assumed to correspond to the luminosity of 
observed ultraluminous X-ray sources, then 



Thus, an IMBH has no time to increase appreciably itsmass in the lifetime of theUniverse 
even if it always moves in the disk. A larger mass growth could be expected for the 
Eddington luminosity of a black hole at rest in a molecular cloud. However, the probability 
of an IMBH being captured into the galactic disk is very low. The probability of an IMBH 
being directly in a cloud of cold gas is even lower, because molecular clouds occupy only a 
disk volume fraction / MC ~ 0.017 (Mapelli et al. 2006). Consequently, no blackhole mass 
redistribution occurs and the region of the IMBH mass function remains suppressed. 



Here, we presented a new SMBH formation mechanism in which the formation of 
IMBHs is suppressed. Whereas the existence of SMBHs at galactic centers is already 




(11) 



CONCLUSIONS 
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beyond doubt, the currently available observational data provide only circumstantial evi- 
dence for the existence of IMBHs. Since in some models for the formation of black holes 
their mass spectrum is close to a power law, the formation of SMBHs is inevitably ac- 
companied by the formation of a large number of IMBHs. In the model we developed, 
this does not happen. On the contrary, the formation of IMBHs is suppressed through 
the classical motion of the phase of a scalar field additional to the inflaton. If detailed 
astronomical observations will establish that the number of IMBHs in the Universe is very 
small, then this fact can be explained in the presented model. 
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Figure 1: Dependence of 9 on the number of e-folds at A = 1.75H and / = 10H: the solid 
line and the dots represent the numerical calculation and the analytical calculation in the 
slow-rolling approximation, respectively; the fitting straight line is indicated by dashes, its 
slope is a = 0.028. 
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Figure 2: PBH mass distribution: the solid line indicates the distribution for a = 0.028 
and A 1.75if; the dashed line indicates the distribution for a = 0. 



